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Abstract 

The  molecular  dynamics  (MD)  computer  simulation  technique  has  been  used  to  study  electrolyte/cathode  interfaces  formed  in 
Li-based  thin  film  oxide  solid  state  ionic  devices  at  the  atomistic  level.  The  solid  electrolytes  are  lithium  silicate  glasses  while  the 
cathodes  are  V205  or  W03  crystals.  The  work  presented  in  this  paper  will  focus  on  the  behavior  at  the  glass/ V,05  interface.  The  MD 
simulation  technique  has  been  successfully  used  to  simulate  a  variety  of  silicate  glasses  and  glass  surfaces,  with  results  consistent  with  a 
variety  of  experimental  data.  The  simulations  of  the  vanadia  crystal  reproduce  the  experimental  crystal  structures,  vibrational  frequency, 
and  the  appropriate  phase  transition  of  V205  as  Li  ions  enter  the  crystal.  The  simulations  have  also  shown  that  Li  transport  into  the  crystal 
is  affected  by  the  orientation  of  the  crystal  at  the  interface  as  well  as  by  surface  roughness.  While  the  crystal  oriented  with  the  (001) 
planes  parallel  to  the  crystal/glass  interface  shows  the  appropriate  phase  transition  to  the  S-LiV,05  phase  as  Li  ions  enter  the  crystal,  the 
work  presented  here  shows  that  the  crystal  oriented  with  the  (100)  planes  parallel  to  the  interface  does  not  transform.  The  difference  is 
attributed  to  the  effect  of  interface  bonding  between  the  ions  in  the  first  crystal  layer  and  those  in  the  glass  surface.  The  simulations  show 
a  relaxation  occurring  in  a  lithium  metasilicate  glass  electrolyte  but  not  in  a  lithium  disilicate  electrolyte.  In  addition,  relaxation  at  the 
interface  between  a  roughened  glass  surface  and  the  crystal  creates  a  distortion  in  the  crystal  planes  in  immediate  contact  with  the  glass 
that  creates  an  induced  strain  in  the  crystal.  ©2000  Elsevier  Science  S.A.  All  rights  reserved. 
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1.  Introduction 

Solid  state  lithium  ion  batteries  have  been  considered 
for  microbattery  applications  for  several  decades  [1-6]  and 
new  materials  are  constantly  being  evaluated  [7-14].  While 
bulk  properties  of  the  individual  components  are  fairly 
well  understood,  the  properties  at  the  different  interfaces 
that  are  created  in  these  systems  are  not  as  well  under¬ 
stood.  The  difficulty  in  obtaining  experimental  data  of 
interface  structure  and  behavior  at  an  atomistic  level  has 
prompted  the  application  of  computational  techniques  to 
study  these  interfaces. 
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The  molecular  dynamics  (MD)  computer  simulation 
technique  is  one  important  computational  approach  to  ad¬ 
dressing  these  problems.  Classical  MD  simulations  rest  on 
the  ability  of  the  assumed  interatomic  potential,  from 
which  the  force  between  atoms  is  derived,  to  be  reasonably 
accurate  in  its  description  of  atomistic  behavior.  In  the 
simulations  performed  here,  the  system  consists  of  lithium 
silicate  electrolyte  glasses  and  vanadia  cathodes. 

We  have  performed  a  large  number  of  simulations  of 
silicate  glasses,  glass  surfaces,  crystals,  and  molecules 
[15-33]. 

Using  a  multibody  potential,  the  simulations  generated 
the  structure  of  the  simulated  silica  glass  similar  to  the 
experimental  structure  function  [18].  The  structure  of  alkali 
silicate  glasses  [24]  gave  the  maximum  in  the  first  peak  in 
the  pair  distribution  function  (PDF)  for  Na-O  at  2.36  A, 
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similar  to  the  extended  X-ray  absorption  fine  structure 
(EXAFS)  data  [34].  More  interestingly,  the  simulations 
showed  the  cause  for  the  need  for  a  two-shell  model  in  the 
curve-fitting  of  the  EXAFS  data.  Thus,  the  simulations 
reproduced  the  overall  EXAFS  curve  for  Na-O  exactly, 
and  showed  which  pairs  caused  the  first  shell  in  the 
EXAFS  calculations. 

Simulations  of  silica  and  alkali  silicate  surfaces  matched 
experimental  data,  without  modifying  the  potentials  to 
specifically  address  surface  properties.  Simulations  of  dry 
[19]  and  wet  silica  [21]  surfaces  using  this  multibody 
potential  resulted  in  structures  that  were  consistent  with 
experimental  data  (not  only  qualitatively,  but  also  quantita¬ 
tively).  Early  simulations  of  alkali  silicate  glass  surfaces 
showed  preferential  relaxation  of  K  and  Na  ions  to  the 
outer  glass  surface,  but  no  such  behavior  for  Li  ions 
[17,35,36],  consistent  with  ion  scattering  spectroscopy  (ISS) 
studies  of  such  glasses  [37].  While  these  early  studies  used 
pair  potentials,  similar  results  are  observed  in  current 
simulations  with  the  multibody  potentials. 

Simulations  of  the  bulk  and  the  surfaces  of  a-alumina 
and  y-alumina  [28]  were  done  using  multibody  potentials 
which  allow  for  the  two  possible  coordinations  of  the  A1 
atoms  in  the  alumina  crystals  (i.e.,  the  A1  ions  can  be 
either  octahedrally  coordinated  in  a-alumina  and  tetrahe- 
drally  and  octahedrally  coordinated  in  y-alumina).  The 
simulation  results  of  the  alumina  systems  are  in  agreement 
with  experimental  X-ray  diffraction  (XRD)  data  and 
molecular  static  calculations.  The  surface  energy  of  the 
Al-terminated  (0001)  alumina  surface  was  found  to  be  2.0 
J/m2,  consistent  with  the  more  recent  DFT  calculations 
which  put  the  value  at  1.95  J/m2  [38].  In  addition,  the 
simulations  show  a  relaxation  of  the  A1  ions  in  the  Al- 
terminated  (0001)  surface  inwards,  similar  to  the  ab  initio 
calculations. 

In  simulations  relevant  to  thin  film  batteries,  we  have 
simulated  both  V-,05  and  W03  crystals  in  contact  with 
lithium  silicate  glasses  using  multibody  potentials  in  the 
MD  technique  [32,39]. 

In  the  more  recent  simulations  of  the  cathode/glass 
interface,  emphasis  has  been  on  the  V,0-  cathode.  Con¬ 
stant  pressure  simulations  of  the  crystal  structure  matched 
XRD  patterns  and  unit  cell  parameters  from  the  literature, 
as  well  as  the  vibrational  spectrum  [39].  Charge  variation 
in  the  V  ions  caused  by  the  local  presence  of  the  Li  was 
incorporated  in  the  simulations  [39]. 

The  simulations  showed  a  phase  transformation  from 
V2Os  to  8-LiV-,05  as  Li  ions  enter  the  crystal,  which 
follows  the  experimental  phase  diagram  [39].  Interestingly, 
the  interatomic  potentials  were  developed  for  V205  and 
y-LiV205,  not  8-LiV205.  However,  the  8-LiV905  phase  is 
stable  with  this  multibody  potential,  providing  an  indica¬ 
tion  of  its  value  in  these  simulations.  The  phase  transfor¬ 
mation  caused  by  the  presence  of  the  Li  ions  in  the 
vanadia  occurs  by  a  glide  of  the  (001)  layers  in  the  (010) 
direction.  This  phase  transformation  was  observed  for  the 


V205  crystal  oriented  with  the  (001)  planes  parallel  to  the 
glass  surface  (called  the  ‘(001)  orientation'),  but  not  when 
the  V205  was  oriented  with  the  (010)  planes  parallel  to  the 
glass  surface  (called  the  ‘(010)  orientation')  [39].  The 
transformation  was  shown  to  be  a  lateral  shift  occurring 
between  the  layers  of  the  cathode  crystal  due  to  the 
presence  of  lithium  between  those  layers  (see  Fig.  1). 
Those  upper  layers  that  did  not  have  Li  ions  between  them 
did  not  show  the  glide  transformation  (Fig.  1),  but  did 
transform  after  additional  Li  ions  migrated  to  those  layers 
[39], 

We  attribute  the  lack  of  transformation  of  the  (010) 
orientation  to  the  required  glide  of  the  (001)  layers  in  the 
(010)  direction,  which  is  physically  inhibited  when  the 
(010)  planes  are  parallel  to  the  surface  (see  Fig.  2). 
However,  with  the  (100)  planes  parallel  to  the  surface, 
glide  in  the  (010)  direction  is  not  inhibited  (see  Fig.  2), 
but  would  be  parallel  to  the  interface.  This  orientation  was 
studied  in  recent  simulations  and  the  results  are  presented 
in  this  paper. 

Results  of  our  simulations  of  the  (lOO)-oriented  crys¬ 
talline  vanadia  in  contact  with  the  electrolyte  glass  in  order 
to  evaluate  the  transformation  to  the  8-LiV905  phase  are 


Fig.  1.  The  (001)  planes  of  vanadia  oriented  parallel  to  crystal/glass 
interface.  View  is  in  the  (100)  direction  and  shows  the  shift  between 
layers  1  and  2  caused  by  the  presence  of  Li  ions  between  the  layers.  No 
shift  observed  in  the  upper  layers. 
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Fig.  2.  Three  orientations  schematically  showing  (010)  glide  directions  necessary  for  the  phase  transformation  to  the  delta  lithium  vanadate  phase. 


reported  below.  In  addition,  Li  buildup  at  the  interface  as  a 
function  of  orientation  is  also  discussed. 


2.  Computational  procedure 

2.1.  MD  technique 


The  details  of  the  computational  procedure  have  been 
previously  presented  [39,40],  but  important  aspects  will  be 
briefly  reiterated  here.  The  classical  MD  simulations  in¬ 
volve  solving  Newton's  equations  of  motion  for  a  system 
of  atoms  (or  ions)  interacting  via  an  assumed  interatomic 
potential.  The  general  form  of  the  interatomic  multibody 
potential  used  here  is  given  as: 

V(r)  =  £  v2(ri,rj)  +  £  f3( rt,rj,rk) . 

i<j  i.j.k 


The  two-body  potential  (V2)  is  given  as  a  modified 
Born-Mayer-Huggins  ionic  potential: 

<h<lje2 
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with  rjj  =  \ri  —  rj\,  q  are  the  charges,  ply-,  Ai;-,  /3(J  are 
adjustable  parameters,  £  is  the  complementary  error  func¬ 
tion  that  acts  to  reduce  the  charges  as  a  function  of 
separation  distance,  and  VASF  is  used  to  add  additional 
structure  to  the  potential  and  is  given  as: 

6 

vr 


=  E 


x= !  l+sMbu(rij~Cu))  ' 


VjJSF  is  zero  for  use  in  the  silicates,  but  was  useful  in  the 
simulations  of  vanadia  and  for  systems  containing  H  ions, 
as  described  below. 


The  three-body  potential  is  given  as: 
v3 (ri’rj’rk)  =  <£„('); -'A  •%;) 

=  *jikexP[Vij/(rij  ~ri)  +  y-k/(  rik  ~  lc)] 

XDyit 

if  ru  <  r,c  or  rik  <  r,c, 

=  <i>jik(rij’rikJjik)  =  0  if  r,j  >  rt  or  rlk  >  r£, 

f3ijk ,  A jjk,  rc ,  and  6J[k  are  adjustable  parameters  and  0jjk 
is  the  angle  between  the  vectors  r-  —  r,  and  rk  —  r,.  The 
three -body  contribution  iljik  comes  from  the  atoms  i,  j, 
and  k,  where  atoms  j  and  k  are  covalently  bonded  to  the 
central  atom  i.  The  strength  of  the  three -body  term  can  be 
varied  so  as  to  reflect  the  covalent  nature  of  the  bonds  in 
question  (or  set  to  zero  if  desired).  Parameters  for  the 
potentials  have  been  previously  presented  [39]. 

A  fifth-order  Nordsieck-Gear  predictor-corrector  algo¬ 
rithm  was  used  to  solve  Newton’s  equation  of  motion. 
Both  NVE  (microcanonical  ensemble)  and  NpT  (canoni¬ 
cal)  simulations  were  performed.  NpT  simulations  were 
run  using  a  modified  Berendsen  algorithm  [41]  to  take  into 
account  non-isotropic  conditions  [28].  Adjustments  to  the 
previously  published  procedure  [39,40]  were  made  only  in 
the  manner  by  which  surface  and  interfaces  were  formed 
and  are  presented  as  needed  below. 

2.2.  Surface  and  interface  formation 

The  interface  between  the  cathode  and  electrolyte  was 
formed  normal  to  the  s-axis  by  placing  a  V205  crystal 
surface  in  contact  with  a  surface  of  a  lithium  metasilicate 
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(LizO  -  Si  02)  glass.  System  sizes  were  around  21-25  A 
in  X  and  Y  (parallel  to  the  surface),  with  a  Z  dimension 
~  25  A  for  the  glass  phase  and  ~  45  A  for  the  vanadia 
phase.  Approximately  2800  ions  were  used  and  multiple 
simulations  were  performed.  The  sizes  varied  because  of 
the  difference  in  the  orientation  of  the  crystals  used  in  the 
study.  The  glass  surface  farthest  from  the  interface  was 
terminated  with  frozen  atoms  and  the  crystal  surface  far¬ 
thest  from  the  interface  was  terminated  with  free  space  in 
order  to  simulate  a  thin  film  of  vanadia  on  an  electrolyte 
substrate. 

Both  the  (001)  and  the  (100)  crystallographic  orienta¬ 
tions  in  contact  with  the  glass  were  simulated.  For  the 
(001)  orientation,  two  surface  terminations  were  created: 
one  surface  was  terminated  with  vanadium  atoms  and  the 
other  surface  was  terminated  with  the  oxygen  atoms  at  the 
end  of  the  vanadyl  bond.  Both  V  and  O  terminations  were 
also  used  for  the  (lOO)-oriented  systems.  Each  termination 
was  made  by  cleaving  the  crystal  normal  to  the  direction 
of  the  appropriate  orientation  at  a  position  which  produced 
the  desired  vanadium  rich  or  oxygen-rich  surface.  The 
newly  cleaved  crystal  surfaces  were  not  annealed.  Fig.  3 
shows  the  crystal  terminations  for  both  the  (00 1)  and  (100) 
orientations. 

The  lithium  metasilicate,  Li20  •  Si02,  glass  was  made 
using  a  constant  volume  and  energy  (NVE)  melt/quench 
cycle.  The  initial  simulation  box  size  matched  the  room 
temperature  X  and  Y  dimensions  of  the  V905  crystal  of 
the  appropriate  orientation,  but  was  scaled  in  accordance 
with  a  thermal  expansion  coefficient  of  120.0  X  10_7oC_1 


at  elevated  temperatures.  The  density  of  the  glass  at  300  K 
was  2.323  g/cm. 

The  periodic  boundary  condition  in  the  Z  dimension  of 
the  simulation  box  was  removed  to  create  free  glass  sur¬ 
faces  perpendicular  to  the  Z  axis.  Finally,  the  glass  system 
was  put  through  an  annealing  cycle  in  order  to  relax  the 
newly  formed  surfaces  of  broken  and  strained  bonds.  The 
X  and  Y  dimensions  for  the  individual  glass  and  crystal 
systems  were  similar  to  the  previous  studies,  but  thicker 
crystals  (larger  Z  dimension)  were  used  here. 

Once  the  glass  and  crystal  systems  had  been  indepen¬ 
dently  formed,  they  were  joined  together  in  one  simulation 
box  to  form  an  ‘Interface  System’.  The  X  and  Y  dimen¬ 
sions  of  the  simulation  box  were  kept  the  same  as  the 
crystal  simulation  box.  In  the  Interface  System,  the  elec¬ 
trolyte  glass  was  annealed  one  more  time  to  account  for 
the  presence  of  the  newly  added  cathode  material  above  its 
top  surface.  During  this  annealing,  the  cathode  ions  were 
kept  immobilized  and  all  glass  ions  were  allowed  to  move. 
The  relatively  short  time  annealing  was  done  sequentially 
at  the  following  temperatures  and  times:  300  K  for  0.5  ps, 
1000  K  for  1  ps,  600  K  for  0.5  ps,  and  300  K  for  0.5  ps. 
The  normal  long-time  simulations  followed  annealing. 

For  the  interface  simulations,  an  imaginary  plane,  called 
the  ‘surface  plane’,  was  set  above  the  electrolyte  surface  to 
mark  the  point  at  which  lithium  ions  were  considered  to 
have  left  the  glass  and  entered  the  interface.  For  each 
lithium  that  crossed  this  surface  plane,  a  new  lithium  was 
added  to  the  bottom  of  the  electrolyte  in  a  volume  above 
the  frozen  atoms.  Adding  lithium  in  this  manner  main- 


Fig.  3.  The  O  and  V  terminations  of  the  (100)  termination  planes. 


194 


S.H.  Garofalini,  P.  Shadwell  /Journal  of  Power  Sources  89  (2000)  190-200 


tained  the  glass  stoichiometry  and  mimicked  a  supply  of 
lithium  from  an  imaginary  anode. 

To  find  the  position  of  the  ‘surface  plane’,  a  density 
profile  of  the  electrolyte  was  generated  along  the  7-axis,  as 
shown  in  Fig.  4.  From  the  plot,  the  point  at  which  the 
density  of  the  electrolyte  falls  to  roughly  one-half  the  bulk 
value  —  the  Gibbs  Dividing  Surface  —  was  chosen  to  be 
the  location  of  the  ‘surface  plane’. 

The  separation  distance  between  the  glass  and  the  crys¬ 
tal  was  determined  by  reference  to  the  distance  between 
the  ‘surface  plane’  and  the  bottom  layer  of  the  crystal. 
Separation  distances  of  2,  3,  and  4  A  were  used  in  order  to 
determine  the  effect  of  this  feature  on  interface  structure 
and  Li  intercalation.  As  a  result  of  this  aspect  of  the  work, 
additional  studies  were  performed  in  order  to  address 
relaxation  of  the  glass  surface  in  the  presence  of  the 
crystal,  as  will  be  presented  below. 

After  the  crystal  and  glass  were  joined  and  annealed, 
simulation  runs  at  300  K  for  varying  times  (from  10  to  100 
ps)  under  NpT  conditions  were  used  in  order  to  gather  the 
statistics  of  interface  behavior.  All  ions  except  those  in  the 
lower  frozen  region  of  the  glass  were  allowed  to  move. 

As  mentioned  above,  studying  the  effect  of  the  initial 
separation  distance  between  the  crystal  and  the  glass  when 
they  were  first  brought  together  led  to  an  additional  study. 
Specifically,  a  large  glass  ( ~  46  X  ~  64  A2  in  X  and  Y ) 
was  put  within  2  A  of  a  three -layer  crystal.  The  center 
third  of  the  glass  in  the  largest  dimension  had  glass  atoms 
at  the  surface  removed  so  that  this  region  was  ~  4-5  A 
away  from  the  crystal.  Both  a  lithium  metasilicate,  Li  20  • 
Si02,  glass  and  a  lithium  disilicate,  Li90  •  2Si02,  glass 
were  used  in  this  study.  The  metasilicate  contains  more  Li 
ions  and  shows  considerably  less  network  connectivity 
(Si — O — Si  bond  connectivity)  than  does  the  disilicate.  The 
effect  of  this  compositional  change  on  smoothing  of  the 
glass  surface  towards  the  crystal  was  evaluated.  It  also 
provided  additional  information  regarding  interface  behav¬ 
ior  that  will  be  presented  below. 


DISTANCE  PERPENDICULAR  TO  INTERFACE  (A) 

Fig.  4.  Density  profile  of  electrolyte  glass.  The  x  marks  values  along 
z-axis,  which  was  used  to  define  the  ‘surface  plane’. 


3.  Results  and  discussion 


3.1.  Transformation 

At  each  separation  distance,  the  systems  with  the  elec¬ 
trolyte  joined  to  the  (001)  surface  of  the  cathode  showed 
that  lithium  intercalation  caused  a  transformation  of  the 
cathode  from  V205  to  8-LiV-,05.  However,  the  termination 
location  of  this  orientation,  O  or  V  termination,  altered  the 
location  of  the  observed  transformation. 

With  the  O  termination  of  the  (00 1)  orientation,  the 
transformation  occurred  between  the  first  and  second  crys¬ 
talline  layers  as  Li  filled  the  sites  between  these  layers, 
similar  to  that  shown  previously  [39,40]  and  shown  in  Fig. 
5.  However,  when  the  cathode  in  this  orientation  was 
terminated  with  V  ions,  the  transformation  was  prohibited 
from  occurring  between  the  first  and  second  layers  of  the 
cathode  and  only  occurred  after  Li  entered  between  the 
second  and  third  crystal  layers,  as  shown  is  Fig.  5.  A  better 
view  of  this  interface  is  shown  in  Fig.  6,  where  the 
V— O— V  bonding  between  the  first  and  second  layers  is 
obvious.  This  bonding  prevents  the  necessary  glide  be¬ 
tween  the  first  and  second  layers.  The  V  whose  vanadyl 
oxygen  bonds  were  removed  in  the  first  layer  are  underco¬ 
ordinated.  Combined  with  the  repulsive  forces  from  the  Li 
ions  in  the  glass  which  quickly  move  into  the  interface, 
these  V  ions  interact  with  the  downward  facing  vanadyls 
from  the  second  layer  to  which  they  bond.  A  shortening  of 
the  long  ‘bond’  from  2.791  A  to  roughly  1.8  A  occurs  with 
a  concomitant  elongation  of  the  vanadyl  bond  from  the 
second  layer  from  1.4  to  1.8  A.  The  similarity  of  these  two 
V— O  bonds  is  shown  in  Fig.  6  and  is  similar  to  a  bond  in 
V203.  Once  the  second  crystal  layer  is  bonded  to  the  first 
crystal  layer,  the  two  are  restricted  from  lateral  movement 
and  hence,  no  transformation  occurs  between  these  layers, 
even  though  Li  ions  are  present. 

While  we  would  expect  the  O-terminated  (001)  plane  to 
be  the  more  normal  situation  structurally,  the  study  with 
the  V-terminated  surface  in  contact  with  the  glass  provides 
for  an  interesting  speculation.  If  some  portion  of  a  V205 
crystal  shows  sufficient  interlayer  bond  formation  (for  any 
reason,  such  as  grain  boundaries  or  point  defects,  or  loss  of 
oxygen  to  locally  reduce  O  concentration,  etc.,  not  just  that 
presented  here  in  the  simulations),  one  could  expect  a 
change  in  transformation  behavior  as  well  as  diffusive 
behavior  of  the  Li  ions.  As  presented  previously  [40],  Li 
diffusion  in  V-,05  is  very  anisotropic.  Simulation  results 
showed  that  the  V205  crystal  has  a  minimum  activation 
barrier  of  0.87  eV  for  Li  migration  in  the  (010)  direction 
and  a  much  higher  barrier  of  2.47  eV  for  diffusion  in  the 
(001)  direction.  For  the  8-LiV-,05  phase,  a  minimum 
barrier  of  0.81  eV  in  the  (010)  direction  was  observed  vs. 
1.79  eV  in  the  (001)  direction  [40].  The  ‘cage’  set  up  by 
the  bonds  shown  in  Fig.  6  would  be  expected  to  slow  one 


O  TERMINATED  V  TERMINATED 


Fig.  5.  Phase  transformation,  labeled  PT  in  the  figure,  occurs  by  sliding  between  layers  1  and  2  in  the  O-terminated  surface,  as  shown  in  Fig.  1.  In  the  V-terminated  surface,  interplanar  bonding,  labeled  IB, 
prevents  the  glide  between  layers  1  and  2  with  Li  present,  but  layers  2  and  3  do  show  the  delta  phase  shift.  The  IB  is  shown  better  in  Fig.  6. 
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Fig.  6.  Slide  view  of  V-terminated  (001)  interface. 


of  the  more  rapid  diffusion  directions  and  alter  Li  trans¬ 
port,  perhaps  even  trapping  some  Li  ions  in  the  crystal 
during  cycling  in  a  manner  otherwise  unanticipated.  This 
heterogeneity  in  activation  barriers  for  Li  diffusion  in 
vanadia  also  has  some  important  implications  with  respect 
to  Li  transport  across  the  cathode/glass  interface,  as  pre¬ 
sented  below. 

Although  the  crystal  with  the  (100)  orientation  could 
glide  parallel  to  the  crystal /glass  interface,  as  indicated  in 
Fig.  3,  the  crystal  with  this  orientation  did  not  transform  to 
the  delta  phase,  regardless  of  the  separation  distance.  This 
lack  of  phase  transformation  is  similar  to  the  results  for  the 
(OlO)-oriented  crystal  [39].  The  absence  of  any  transforma¬ 
tion  in  the  (lOO)-oriented  system  is  attributed  to  interface 
bonding  between  the  crystal  and  the  glass.  Within  20  ps, 
approximately  20-50%  of  the  crystal  ions  at  the  interface 
bonded  to  the  network  ions  in  the  glass.  The  higher  percent 
of  bonding  occurred  at  the  smaller  initial  separation  dis¬ 
tance,  while  the  lower  bonding  percent  occurred  at  the 
larger  separation  distance.  Such  bonding  inhibits  the  slid¬ 
ing  of  the  individual  crystal  layers.  So  while  the  crystalline 
layers  in  the  (100)  orientation  could  slide  over  the  elec¬ 
trolyte  surface,  as  shown  schematically  in  Fig.  3,  the 
bonding  between  the  crystal  and  the  glass  prevents  this  and 
no  transformation  is  observed. 

Such  bond  formation  in  the  (001)  orientation  had  no 
affect  on  the  upper  (001)  planes.  Thus,  even  though  bond 
formation  could  occur  between  the  first  layer  of  ions  in  the 
crystal  and  the  glass  surface  in  the  (OOl)-oriented  system, 
thus  anchoring  this  first  layer  to  the  glass,  the  upper 
crystalline  layers  were  unaffected  by  this  bonding.  There¬ 
fore,  when  sufficient  Li  ions  entered  between  layers  of  the 
(OOl)-oriented  crystal,  the  upper  layer  could  glide  in  the 
(010)  direction,  affecting  the  transformation  to  the  8- 
LiV205  phase. 


3.2.  Interface  relaxation 

While  the  different  separation  distances  between  the 
glass  and  the  crystal  at  the  start  of  the  simulation  did  not 
have  a  major  effect  on  the  observed  phase  transformation 
(or  lack  thereof  in  the  (100)  orientations),  there  was  an¬ 
other  interesting  result  observed  in  the  simulations.  The 
crystal  moved  slightly  inward  or  outward,  depending  on 
orientation  and  termination.  This  was  anticipated  under 
these  constant  pressure  simulations  which  had  a 
crystal/vacuum  interface  away  from  the  glass  interface 
that  enabled  the  crystal  to  ‘float’  above  the  glass  and  reach 
some  equilibrium  spacing  with  the  glass.  Of  course,  this 
would  be  affected  by  the  interactions  at  the  crystal/glass 
interface,  such  as  bonding  and  Li  migration. 

More  importantly,  in  addition  to  the  obvious  motion  of 
Li  ions  towards  the  crystal,  the  simulations  also  showed 
movement  of  the  network  ions  (Si  and  O)  in  the  lithium 
metasilicate  glass  towards  the  crystalline  layer  in  the  sys¬ 
tems  with  the  larger  initial  separation  distance.  Since  the 
system  is  periodic  in  the  X  and  Y  directions  (parallel  to 
the  interface),  movement  of  the  glass  towards  the  crystal 
could  occur  with  no  restraint  laterally.  Since  the  glass  was 
set  at  a  particular  initial  separation  distance,  the  whole 
glass  surface  would  experience  relatively  similar  forces  at 
the  interface.  Flowever,  similar  simulations  using  lithium 
disilicate  and  lithium  trisilicate  glasses  showed  no  such 
relaxation  of  the  network  ions  in  these  glasses  towards  the 
crystal.  The  disilicate  and  trisilicate  glasses  contain  far 
fewer  Li  ions  than  in  the  metasilicate  glass,  thus  giving 
glasses  with  more  network  connectivity.  The  greater 
three-dimensionality  of  the  Si— O  connectivity  present  in 
the  di-  and  trisilicate  glasses  creates  a  less  mobile  network. 

These  results  raise  the  interesting  question  regarding  the 
response  of  the  glassy  phase  if  the  glass  surface  were  not 
uniformly  separated  from  the  crystal.  That  is,  what  would 
be  the  response  of  the  glassy  phase  if  the  surface  had 
regions  with  different  separation  distances  from  the  crys¬ 
tal?  The  results  would  relate  to  cases  where  the  glass 
surface  is  topographically  much  rougher  than  the  normal 
atomistic  roughness  inherent  for  silica  (silicate)  glass  sur¬ 
faces  [23]. 

In  order  to  address  this  question,  additional  simulations 
were  performed  on  very  large  systems  containing  nearly 
10,000  ions.  A  lithium  metasilicate  glass  was  created  with 
a  Y  dimension  three  times  larger  than  the  previous  work, 
an  X  dimension  two  times  larger  than  the  previous  work. 
The  crystal/glass  separation  distance  was  set  at  2.0  A. 
However,  an  additional  2  A  of  glass  ions  in  the  middle 
third  of  the  glass  in  the  Y  dimension  was  removed  from 
the  glass  surface,  creating  a  channel  in  the  glass  surface,  as 
shown  in  the  schematic  below.  This  created  a  section  of 
the  glass  that  had  a  4.0-4. 5  A  separation  distance  between 
the  glass  and  the  crystal  in  this  portion  of  the  interface.  A 
schematic  of  the  surface  of  this  glass  is  shown  in  Fig.  7. 
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Fig.  7.  Schematic  drawing  of  the  glass  surface  with  a  channel  of  glass 
ions  removed  from  the  middle  third  in  the  Y  dimension. 


Fig.  8a  shows  a  snapshot  of  the  configuration  looking 
from  the  (100)  direction  (the  horizontal  dimension  is  Y). 
The  crystal  is  in  the  (001)  orientation,  with  a  V  termination 
plane.  The  sizes  of  the  ions  are  only  for  the  graphics,  with 
Li  ions  being  the  largest  spheres  in  the  figure.  Si— O  and 
V— O  bonds  are  drawn  in  the  figure.  Note  the  larger 


separation  distance  in  the  middle  third  of  the  image  in 
comparison  to  the  outer  third,  similar  to  the  schematic 
drawing  in  Fig.  7.  Because  of  the  inherent  atomic  rough¬ 
ness  of  the  silicate  surface,  the  actual  atomistic  structure 
shown  in  Fig.  8  is  not  as  clean  (straight  lines)  as  the 
schematic  in  Fig.  7.  Fig.  8b  shows  the  interface  after  10  ps 
at  300  K.  The  ions  in  the  middle  third  of  the  glass  have 
moved  upward  towards  the  crystal  and  smoothened  the 
interface  considerably  in  comparison  to  the  initial  configu¬ 
ration. 

A  similar  simulation  using  the  lithium  disilicate  (LizO  • 
2SiCL)  glass  resulted  in  little  motion  of  the  glass  towards 
filling  the  crystal/glass  gap.  The  higher  network  connec¬ 
tivity  of  the  disilicate  vs.  the  metasilicate  inhibited  the 
relaxation  of  the  former  in  comparison  to  the  relaxation 


Fig.  8.  Lithium  metasilicate  with  three  layers  of  vanadia  crystal  on  top.  (a)  Starting  configuration  showing  2  A  separation  between  the  crystal  and  the  glass 
in  the  1st  and  3rd  of  the  glass  (as  measured  in  the  horizontal  direction),  with  additional  glass  ions  removed  from  the  glass  surface  in  the  middle  third  of  the 
glass. 
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LITHIUM  METASILICATE  SURFACE 


Y  Position  (A) 

LITHIUM  DISILICATE  SURFACE 

36.5 


Y  Position  (A) 

Fig.  9.  Relaxation  of  ions  in  glass  towards  crystal  in  ‘channel’  structure 
(see  Fig.  7).  Maximum  Z  coordinate  of  glass  ions  averaged  over  X  for 
each  Y  value,  (a)  Lithium  metasilicate,  (b)  lithium  disilicate.  (Note 
difference  in  ordinate  between  a  and  b.) 


that  was  seen  in  the  metasilicate  case.  The  simulation  with 
the  disilicate  glass  was  run  for  over  100  ps  with  no 
additional  relaxation  of  the  glass. 


Additional  representations  of  the  glass  surface  relax¬ 
ation  behavior  are  shown  in  Fig.  9.  Here,  the  highest  Z 
positions  of  atoms  in  the  glass  surface  as  a  function  of  the 
atom’s  Y  position,  averaged  over  all  X  positions  at  that  Y 
value,  are  plotted  for  the  initial  glass  and  the  final  configu¬ 
ration.  Clearly,  in  Fig.  9a,  the  middle  third  of  the  metasili¬ 
cate  glass  relaxes  upward  (towards  the  crystal  which  is  not 
shown),  while  in  Fig.  9b,  the  middle  third  of  the  disilicate 
glass  shows  little  such  motion.  Because  of  the  greater 
amount  of  network  connectivity  in  the  disilicate  glass  than 
the  metasilicate  glass,  there  is  a  greater  restriction  in  the 
motion  of  the  silica  chains  and  rings. 

Since  the  outer  thirds  (first  third  and  third  third)  of  the 
metasilicate  glass  is  only  2  A  away  from  the  crystal  and 
shows  little  motion  towards  the  crystal  (see  Fig.  9a),  the 
connectivity  of  the  middle  third  to  the  two  outer  thirds 
might  be  expected  to  retard  relaxation  of  the  middle  third 
towards  the  crystal.  However,  the  significant  amount  of  Li 
ions  in  the  metasilicate  glass  disrupts  the  network  structure 
so  much  that  short  chains  and  non-bridging  oxygen  domi¬ 
nate  the  structure,  thus  reducing  network  connectivity. 
This  enables  the  middle  third  of  the  glass  to  move  toward 
the  crystal  without  the  restrictive  attachments  to  the  outer 
thirds. 

Also  apparent  in  Fig.  8  of  the  metasilicate  glass  is  the 
buildup  of  Li  ions  at  the  crystal/glass  interface  with  this 
(001)  orientation.  This  buildup  was  not  observed  in  the 
simulations  with  the  (010)  or  (100)  orientations.  Fig.  10 
shows  a  comparison  of  this  difference  in  Li  behavior  at  the 
interface  for  the  (001)-  and  (OlO)-oriented  interfaces.  The 


Fig.  10.  Li  buildup  at  (OOl)-oriented  interface  vs.  no  buildup  in  (010)  orientation. 
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buildup  of  Li  ions  at  the  (001)  interface  can  be  attributed 
to  the  difference  in  the  activation  barriers  to  Li  migration 
in  the  (001)  direction  vs.  that  in  the  (010)  direction.  As 
mentioned  above,  the  barrier  in  the  (001)  direction  is 
~  2.5  eV,  while  the  barrier  in  the  (010)  direction  is  only 
~  0.87  eV  [40].  Activation  barriers  in  8-LiV-,05  are  simi¬ 
larly  anisotropic  [40].  The  activation  energy  for  Li  diffu¬ 
sion  in  the  lithium  silicate  glasses  ranges  from  0.78  to  0.85 
eV  for  various  lithium  silicate  glasses.  Thus,  Li  ions 
diffuse  fairly  rapidly  in  the  glass  towards  the  vanadia 
crystal.  If  the  crystal/glass  interface  has  the  (00 1)  orienta¬ 
tion,  there  is  a  much  higher  activation  barrier  for  Li 
transport  into  the  crystal  in  comparison  to  the  barrier  in  the 
glass,  thus  creating  a  retardation  of  motion  and  a  buildup 
of  Li  ions  at  the  interface.  Since  the  activation  barrier  for 
Li  transport  in  the  (010)  direction  is  not  much  different 
than  that  for  the  glass,  Li  does  not  build  up  at  the 
(OlO)-oriented  interface. 

Analysis  of  the  interface  indicates  that  there  is  a  Li/V 
ratio  of  approximately  1  at  the  interface  in  the  (001) 
orientation.  Experimental  data  in  the  literature  have  been 
interpreted  by  the  authors  to  indicate  a  buildup  of  Li  ions 
at  the  interface  region  [42].  Their  result  of  2.82  (X  10-4) 
C/cm2  would  put  the  Li  concentration  at  ~  1.8  X  10 15 
Li/cm2.  Based  on  the  (001)  orientation,  with  planar  unit 


cell  dimensions  of  11.51  X  3.56  A  ,  a  concentration  of  1 
Li/2  V  in  LiV205  would  create  a  concentration  of  0.98  X 
1015  Li/cm2.  Thus,  their  number  of  ~  1.8  X  1015  Li/cm2 
is  nearly  equivalent  to  1  Li/V  (1.96  X  1015  Li/cm2), 
indicating  a  multilayer  adsorption  of  Li  at  the  interface. 
Our  simulations  show  Li  ions  at  the  interface  equivalent  to 
1  Li/V,  similar  to  the  experimental  data.  In  addition,  as 
mentioned  previously  [39],  the  simulations  showed  Li  con¬ 
centrations  in  the  crystal  layers  with  a  Li/V  ratio  near  1. 

The  simulation  results  regarding  activation  barriers  for 
Li  transport  indicate  that  if  vanadia  crystals  orient  differ¬ 
ently  with  respect  to  the  glass  substrate,  there  will  proba¬ 
bly  be  regions  with  different  conductivity  behavior. 

Finally,  Fig.  8  shows  an  interesting  feature  regarding 
the  structure  of  the  thin  three-layer  crystal  used  in  this 
particular  simulation.  While  in  Fig.  8a,  the  crystalline 
planes  are  straight,  there  is  a  clear  curvature  of  the  planes 
in  Fig.  8b.  This  effect  is  shown  more  clearly  in  Fig.  11, 
which  is  the  same  structure  as  in  Fig.  8,  but  only  a  thin 
section  (into  the  plane  of  the  figure)  is  drawn.  The  crystal 
layers  relax  to  the  shape  of  the  glassy  surface  in  Figs,  lib 
and  8b,  thus  creating  distorted  crystalline  planes  near  the 
interface.  These  crystal  planes  contain  strained  bonds  in 
the  most  distorted  parts  of  the  plane  and  could  affect  the 
mechanical  integrity  of  the  film. 


(a) 


Fig.  11.  Similar  to  Fig.  8,  but  only  a  thin  section  into  the  plane  of  the  figure  shown  here,  (a)  Start  of  run,  (b)  end  of  run,  showing  curvature  of  crystal 
planes. 
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4.  Conclusions 

The  application  of  the  MD  computer  simulation  tech¬ 
nique  to  studies  of  thin  film  lithium  ion  batteries  has 
shown  results  which  are  consistent  with  the  structural 
features  of  the  cathode  crystal  and  the  electrolyte  glass 
films.  The  phase  transformation  to  the  8-LiV,o5  phase  as 
Li  ions  enter  the  crystal  is  affected  by  interfacial  bonding 
between  the  crystal  and  the  glass  in  the  (100)  orientation, 
but  not  the  (001)  orientation  (except  at  the  first  crystal 
layer).  The  large  anisotropy  of  Li  diffusion  into  the  crystal 
alters  migration  of  Li  ions  into  the  crystal,  creating  a 
buildup  of  Li  at  the  interface  in  the  (001)  orientation,  but 
no  such  buildup  in  the  (010)  and  (100)  orientations. 

Relaxation  of  the  ions  in  the  metasilicate  glass  enables  a 
smoothening  of  rough  interfaces,  which  does  not  occur  as 
readily  in  the  lithium  disilicate  glass. 
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